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Abstract 

In the framework of potential NRQCD, we obtain the next-to-leading-log renormalization-group 
running of the matching coefficients for the heavy quarkonium production currents near threshold. 
This allows to obtain S-wave heavy-quarkonium production/annihilation observables with next-to- 
leading-log accuracy within perturbative QCD. In particular, we give expressions for the decays of 
heavy quarkonium to e'^e" and to two photons. We also compute the 0{ma^\T? a) corrections to 
the Hydrogen spectrum. 
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Heavy quark-antiquark systems near threshold are characterized by the small relative ve- 
locity V of the heavy quarks in their center of mass frame. This small parameter produces 
a hierarchy of widely separated scales: m (hard), mv (soft), mv"^ (ultrasoft), ... . The fac- 
torization between them is efficiently achieved by using effective field theories, where one 
can organize the calculation as various perturbative expansions on the ratio of the different 
scales effectively producing an expansion in v. The terms in these series get multiplied by 
parametrically large logs: Inf, which can also be understood as the ratio of the different 
scales appearing in the physical system. Again, effective field theories are very efficient in 
the resummation of these large logs once a renormalization group (RG) analysis of them has 
been performed. This will be the aim of this paper for annihilation and production processes 
near threshold. 

We will restrict ourselves, in this paper, to the situation where Aqcd ^ mal (to be implicit 
in what follows), which is likely to be relevant, at least, for t-i production near threshold. 

NRQCD has an ultraviolet (UV) cut-off vnr = {^p, ^s} satisfying mv <C i'nb. ^ "i. At 
this stage Up ~ Ug. Vp is the UV cut-off of the relative three-momentum of the heavy quark 
and antiquark, p. is the UV cut-off of the three-momentum of the gluons and light quarks. 

Potential NRQCD (pNRQCD) (see [|,|] for details) is defined by its particle content and 
cut-off VpMR = {vpi Vus}-, where Vp is the cut-off of the relative three-momentum of the heavy 
quarks and Vus is the cut-off of the three-momentum of the gluons and light quarks. They 
satisfy the following inequalities: |p| <C i/p -C m and p^/m <^ Vus ^ |p|- Note that no gluons 
or light quarks with momentum of 0(|p|) are kept dynamical in pNRQCD. The motivation 
to integrate out these degrees of freedom is that they do not appear as physical (on-shell) 
states near threshold. Nevertheless, they can appear off-shell and, since their momentum is of 
the order of the relative three-momentum of the heavy quarks, integrating them out produces 
non-local terms (potentials) in three-momentum space. Indeed, these potentials encode the 
non-analytical behavior in the transfer momentum of the heavy quark, k = p — p', of the 
order of the relative three-momentum of the heavy quarks. 

The matching process, which basically means the computation of the potentials, is carried 
out for a given external incoming (outcoming) momentum p (p'). Therefore, one has to sum 
over all of them in the pNRQCD Lagrangian, since they are still physical degrees of freedom 
as far as their momentum is below Vp. In position space, this means that an integral over x, 
the relative distance between the heavy quarks, appears in the Lagrangian when written in 
terms of the heavy quark-antiquark bilinear. 

Within pNRQCD, integrals over p (or x) appear when solving the Schrodinger equation that 
dictates the dynamics of the heavy quarkonium near threshold. At lower orders, these inte- 
grals are finite effectively replacing p by ~ mcts. Nevertheless, at higher orders in quantum 
mechanics perturbation theory and/or if some singular enough operators are introduced (as 
it will be the case of the heavy quarkonium production currents) singularities proportional 
to In z/p appear. These must be absorbed by the potentials or by the matching coefficients of 
the currents. We will describe how to resum the logarithms associated to this cutoff within 
pNRQCD. 
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A RG analysis for non-relativistic systems have been addressed before in Refs. 0-0, where 
they match to an effective theory called vNRQCD. On physical terms, this theory should be 
equivalent to the previously defined pNRQCD once the RG evolution has been performed 
and the soft degrees of freedom have been integrated out, as only ultrasoft gluons and light 
fermions and potential quarks are left as dynamical degrees of freedom. We will compare 
with their results. In some cases disagreement will be found. 

Let us now describe the matching between QCD and pNRQCD within an RG framework. 
For the case where no divergences proportional to In Up appear, the procedure reduces to the 
results of Ref. ^ to which we refer for the notation and background material necessary to 
follow this paper. 

We first address the procedure that gives the running of the potentials. One first does the 
matching from QCD to NRQCD. The latter depends on some matching coefficients: c{i>s) 
and d{vp,Vs)i which can be obtained order by order in as (with Vp = Us) following the 
procedure described in Ref. 0. The c(z/s) stand for the coefficients of the operators that 
already exist in the theory with only one heavy quark (ie. HQET) and the d{i>p, Ug) stand for 
the coefficients of the four heavy fermion operators. The starting point of the renormalization 
group equation can be obtained from these calculations by setting Up = Us = m (up to a 
constant of order one). In principle, we should now compute the running of z/p and z/g. 
The running of the c(z/s) can be obtained using HQET techniques The running of 
the d{i>p, Us) is more complicated. At one-loop, Up does not appear and we effectively have 
d{vp,Vs) ~ divg), whose running can also be obtained using HQET-like techniques |p. At 
higher orders, the dependence on Vp appears and the running of the d{vp, Ug) becomes more 
complicated. Fortunately, we need not compute the running of d in this more general case 
because, as we will see, the relevant running of the d for near threshold observables can be 
obtained within pNRQCD. 

The next step is the matching from NRQCD to pNRQCD. The latter depends on some match- 
ing coefficients (potentials). They typically have the following structure: V{c{h's), d^Up, Ug), ^s, 
v^si r). After matching, any dependence on Ug disappears since the potentials have to be inde- 
pendent of Ug. Therefore, they could be formally written as V{c{l/r), d^Up, 1/r), 1/r, u^g, r). 
These potentials can be obtained order by order in as following the procedure of Refs. pj-Q. 
The integrals in the matching calculation would depend on a factorization scale /i, which 
should correspond either to Ug or to u^g. In the explicit calculation, they could be distin- 
guished by knowing the UV and infrared (IR) behavior of the diagrams: UV divergences are 
proportional to InUg, which should be such as to cancel the Ug scale dependence inherited 
from the NRQCD matching coefficients, and IR divergences to Uug. In practice, however, as 
far as we only want to perform a matching calculation at some given scale fi = Ug = Uug, it is 
not necessary to distinguish between Ug and u^g (or if working order by order in as without 
attempting any log resummation) . 

Before going into the rigorous procedure to obtain the RG equations of the potentials, let 
us first discuss their structure on physical grounds. As we have mentioned, the potential is 
independent of Ug. The independence of the potential with respect Ug allows us to fix the latter 
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to 1/r that, in a way, could be understood as the matching scale for z/^S. Therefore, 1/r, 
the point where the multipole expansion starts, would also provide with the starting point 
of the renormalization group evolution of u^g (up to a constant of order one). The running of 



Uus can then be obtained following the procedure described in Refs. [|11|,§] • At the end of the 
day, we would have V{c{l/r), d^Up, 1/r), 1/r, u^s, r), where the running on Uus is known and 
also the running in 1/r if the d is i/p- independent. So far, the only explicit dependence of the 
potential on Up appears in the d. Nevertheless, the potential is also implicitly dependent on 
the three-momentum of the heavy quarks through the requirement 1/r ~ p -C i/p, and also 
through Uus, since z/^^ needs to fulfill p'^/m -C I'us -C |p|. This latter requirement holds if we 
fix = Vp/m (this constraint tells you how much you can run down in the potential 
before finding the cutoff Vp/m caused by the cutoff of p). 

Within pNRQCD, the potentials should be introduced in the Schrodinger equation. This 
means that integrals over the relative three-momentum of the heavy quarks take place. 
When these integrals are finite one has p ~ 1/r ~ ma^ and p'^/m ^ mal- Therefore, one 
can lower Vus down to ~ mal reproducing the results obtained in Ref. [Q. In some cases, 
in particular in heavy quarkonium creation, the integrals over p are divergent, and the log 
structure is dictated by the ultraviolet behavior of p and 1/r. This means that we can not 
replace 1/r and Vus by their physical expectation values but rather by their cutoffs within the 
integral over p. Therefore, for the RG equation of z/p, the anomalous dimensions will depend 
(at leading order) on V"(c(z/p), d{i'p, Up), z/p, Vp/m, Vp)a and the running will go from Up ^ m 
down to z/p ~ mas- Note that, at this stage, a single cutoff, z/p, exists and the correlation 
of cutoffs can be seen. The importance of the idea that the cutoffs of the non-relativistic 
effective theory should be correlated was first realized by Luke, Manohar and Rothstein in 
Ref. (for an application to QED see |12[)- Note also that at the matching scale z/p ~ m, 
what it would be the ultrasoft cutoff is also of order m. In this sense it should be understood 
the statement in Ref. that ultrasoft gluons appear at the scale m, a point that becomes 
relevant within a RG approach. 

With the above discussion in mind, the matching between NRQCD and pNRQCD could be 
thought as follows. One does the matching by computing the potentials order by order in at 



^ In practice, the potential is often first obtained in momentum space so that one could then set 
= k. Note, however, that this is not equivalent to fix = 1/r, since finite pieces will appear 
after performing the Fourier transform. 
^ Roughly speaking, this result can be thought as expanding In r around In Up in the potential ie. 



V{c{l/r),d{up, 1/r), 1/r, Vp/m, r) ~ V{c{vp),d{vp, i'p),i'p, Vp/m, Vp 

dr 



+ \n{vpr)r—V 



+ ■••■ (1) 

l/T = Vp 



The ln(z/pr) terms may give subleading contributions to the anomalous dimension when introduced 
in divergent integrals over p. The discussion at this stage is not very rigorous and a more precise 
discussion would require a full detailed study within dimensional regularization, which goes beyond 
the aim of this work. Nevertheless, we do not expect it to change the underlying idea, although it 
deserves further investigations. 
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the matching scale i'p = = ^us following the procedure of Refs. (by doing the matching 
at a generic Up some of the running is trivially obtained). The structure of the potential at this 
stage then reads K(c(t'p), (i(z/p, z/p), i/p, i/p, z/p) (and similarly for the derivatives with respect 
Inr of the potential). This provides the starting point of the renormalization group evolution 
of Uus (up to a constant of order one). The running of Uus can then be obtained following 
the procedure described in Refs. ||TT|,|8[|. For the final point of the evolution of u^g, we choose 
Uus = i^p/m. At the end of the day, we obtain \/(c(z/p), (i(z/p, z/p), z/p, I'p/m, Vp) = V{vp). 

The running of Vp goes from Vp = m (this was fixed when the matching between QCD and 
NRQCD was done) up to the physical scale of the problem Vp ~ mag. If the running of the 
NRQCD matching coefficients is known, the above result gives the complete running of the 
potentials. The procedure to get the running of the c is known at any finite order. For the d 
it is just known at one-loop order, since, at this order, it is only z/^- dependent. Nevertheless, 
at higher orders, dependence on Vp appears. Therefore, the above method is not complete 
unless an equation for the running of Vp is provided. This is naturally given within pNRQCD. 
It appears through the iteration of potentials. Let us consider this situation more in detail. 
We first remind what the Hamiltonian in pNRQCD for the singlet sector is (see Ref. |^ for 
notation and further details): 



C4- — - - Uf — - — 



where Cf = {N^ — 1)/ C^Nc) and we will set = C4 = 1 (we will only eventually use C4 for 
tracking of the contribution due to this term). The propagator of the singlet is (formally) 

(3) 



E-h 



At leading order (within an strict expansion in as) the propagator of the singlet reads 



E - /if E - pVm - Cfajr ' 



If we were interested in computing the spectrum at 0{ma^), one should consider the iteration 
of subleading potentials (Shg) in the propagator as follows: 

G,{E)6hsGc{E)---6hsG,{E). (4) 

In general, if these potentials are singular enough, these contributions will produce logarith- 
mic divergences due to potential loops. These divergences can be absorbed in the matching 
coefficients, D^^] and D^g2 g, of the local potentials (those proportional to the ^'•^^(r)) pro- 
viding with the renormalization group equations of these matching coefficients in terms of 
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Vp. Let us explain how it works in detail. Since the singular behavior of the potential loops 
appears for |p| ^ ol^/t, a perturbative expansion in is licit in Gc{E), which can be 
approximated by 

^=^^ =Gf\E)= \ ■ 

— p /m 

Therefore, with the accuracy we aim at in this paper, a practical simplification follows from 
the fact that the Coulomb potential, — C/^, can be considered to be a perturbation as far 
as the computation of the Inz/p ultraviolet divergences is concerned. This means that the 
computation of the anomalous dimension can be organized within an expansion in and 
using the free propagators Gf \ Moreover, each Gf'' produces a potential loop and one extra 
power of m in the numerator, which kills the powers in 1/m of the different potentials. This 
allows the mixing of potentials with different powers in 1/m. One typical example would be 
the diagram in Fig. 1.1 The computation of this diagram would go as follows: 

1 -nCfDf,,^,^^ 1 ^ay^ 1 ttC^/^S 



Using (5(3) (r) = |r = 0)(r = 0|, we can see that the relevant computation reads (instead of 
ay^ one could use as since the non-trivial running of ay^ is a subleading effect. Nevertheless, 
we keep ay^ since it allows to keep track of the contributions due to the Coulomb potentials) 



h=^ 9"! — 1=, 9"! — h" = 0) 

E-p /m r E-p /m 

dV f d'^p m Anay^ m rri^ay^ 1 



{2TiY J {2'KYp'^-mE q2 p2-mE ^ IQn e' 

where D = A + 2e and q = p — p'. This divergence is absorbed in contributing to its 
running at next-to-leading-log (NLL) order as follows 

^v4-Dt^^v) - «v.(^p)^£'(^p) + ■ ■ ■ • (7) 

Therefore, even without knowing the running of the d (which need to be known at NLL 
order in this case), we can obtain the running of the potential (one can also think of trading 
Eq. (1^) into an equation for rf, which is the only unknown parameter within the potential). 
This is so because D)i l is only needed with LL accuracy in the right-hand side of Eq. (0). 



The above method deals with the resummation of logs due to the hard, soft and ultrasoft 
scales. Nevertheless, for some specific kinematical situations even smaller scales could appear. 
Their study, however, goes beyond the aim of this paper. In any case, pNRQCD can be 
considered to be the right starting point to study these kinematical situations. 

3 The diagram in Fig. 1 is also the relevant one in order to obtain the 0{ma^ Ina^) contribution 
to the Hydrogen spectrum. See the Appendix. 
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The matching scale between QCD and NRQCD is t'p ~ ~ m. On the other hand, the 
matching scale between NRQCD and pNRQCD is also the hard scale: z/p ~ z/^/m ~ m. 
Therefore, one could wonder about the necessity of using the intermediate theory NRQCD. 
This is indeed the attitude in Refs. PJT^,|6|J^ , where they directly perform the matching 
between QCD to an effective field theory: vNRQCD that, once the RG evolution has been 
performed and the soft degrees of freedom have been integrated out, should be physically 
equivalent to pNRQCD with Up ~ mas. One motivation for going through NRQCD is that 
it allows to perform the factorization of the hard scale within an effective field theory frame- 
work. In fact, a full factorization of the different regions of momentum that ought to be 
integrated out is achieved within pNRQCD. This extremely simplifies the matching process 
since one deals with only one scale (region of momentum) in the loops at each step. In the 
matching between QCD and NRQCD only hard loops need to be considered, whereas in the 
matching between NRQCD and pNRQCD only soft loops need to be considered. Moreover, 
the structure of the UV cutoffs of the theory is better understood in this way. For instance, 
one can see that all the explicit dependence of the potentials on Up is inherited from the 
d matching coefficients. Within a diagrammatic approach the factorization of the different 



regions of momentum have been achieved using the threshold expansion 113 1. 



Let us now consider the case of the electromagnetic current, which will provide an example 
where to apply the above discussion. The procedure is analogous to the potentials. We first 
do the matching from QCD to NRQCD: 



=&i,nrV^V'x(0) + 0{l/m] 

QCD 



NRQCD 



We will just concentrate in the coefficient &i,nr- Within NRQCD, it should be understood 
as a function of Up and z/^, ie. &i,nr(^'p; i^s)- One should first obtain the matching conditions 
at the hard scale. This has been computed up to two loops [|I^ but we will only need the 



one- loop expression [15|: 



6i,nrKH = 1-2C7^^^, (9) 

TT 

since we only aim to a NLL resummation in this paper. If we compare with the previous 
discussion of the potentials, the matching coefficients d play the role of bi. Therefore, within 
pNRQCD, we will need &i,nr(^'p5 ^'p) = ^i,NR(t'p)- We first have to consider 6i,NR(i^p5 '^s ). In 
this case, unlike for the rf's, there is no running due to Ug at the order of interest. This can 
be easily seen in the Coulomb gauge. Moreover, the matching from NRQCD to pNRQCD 
creates the potentials but let bi unchanged since soft loops or HQET-like calculations give 
zero correction to bi at the order of interest. Formally, 



&i,nrV^V^x(0) 



= 5i,pNR^t^^^(0) 

NRQCD 



(10) 

pNRQCD 



or, in other words, the matching condition reads -Bi,pnr(&i,nr(^'p), '^ms = ^p) = ^i.nrI^'p)- 
The running of u^s is also trivial as there is none at the order of interest (this has to do 
with the fact that we are dealing with an annihilation process). Therefore, we finally have 
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Bi{up) = Bip^^{bi jsf^{h'p),h'p/m) = &i,nr(^'p)- We can see that we are in the analogous 
situation to the running of D^^](z/p) versus the running of d{i'p, Up). We now need the RG 
equation for Bi{vp). This demands to obtain the ultraviolet corrections to the current within 
pNRQCD keeping track of the contributions due to the different potentials. Fortunately, 
this calculation has already been done and we can extract the relevant information from 
Ref. [|l^. The computation goes along the same lines than in the example of Fig. 1. The 
explicit diagrams to be computed for the RG running of Bs{i'p) are given in Fig. 2 (where s 
denotes the spin). From this figure, we can clearly illustrate the structure of the computation. 
0(l/m) corrections to h^^^ only need one potential loop to kill the 1/m coefficient. 0{\/m?) 
corrections to hf^ need two potential loops to kill the coefficient and so on. In the 

situation with more than one potential loop, the additional potential loops can be produced 
without additional 1/m factors coming from the potential only if Coulomb potentials are 
introduced. This explains why the 1/m potential needs zero Coulomb potential insertions, 
the 1/m^ potentials need one Coulomb potential insertions and the 1/m^ term needs two 
Coulomb potential insertions (for the running of D)i 1 and -0^2^, we expect a similar structure). 
In principle, this would be be a never ending story unless there is an small parameter that 
tells us how far we have to go in the calculation in order to achieve some given accuracy. 
This is indeed so. The 1/m potential is a NLL effect p| and therefore higher powers in Z)^^^ 
produce NNLL effects or beyond. On the other hand, the introduction of Coulomb potentials 
brings powers in Os, which suppresses the order of the calculation. In our case, for a NLL 
calculation, the maximum power of the anomalous dimension should be a^. This means that 
with zero insertions (0(l/m) potentials) the matching coefficient {D^p) has to be known 
with NLL accuracy, with one ay^ insertion (0(l/m^) potentials) the matching coefficients 
(Z)(^)) have to be known with LL accuracy and with two cty^ insertions (0(l/m^) potentials) 
the matching coefficients must have no running (this explains why only C4 is considered at 
this order). 



From the above discussion, the RG equation reads @ 



^p-; — Bg = Bs 
a Up 



:ii) 



where Ca = N^, and the RG-improved matching coefficients of the potentials can be read 
from Ref. pl],|8| with the assignment 1/r Vp and Vus Vp/m (see also [[17|,3 for calculations 
of the potentials at finite orders in as) . We have kept the spin explicit so our results will also 
be valid for the pseudoscalar current: 



Q7°75Q(0) 



=feo,NRV^^x(0) + 0(l/m) 

QCD 



(12) 

NRQCD 



^The RG equations of Bg within vNRQCD has been computed in Ref. [^,0]. In principle, they 
are different. Nevertheless, it may happen that fields redefinitions of the potentials may make 
them equal. We have not checked that but it is plausible since, as we will see, for the leading and 
subleading logs (but not beyond) both calculations will agree. 
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with the matching condition [IS]: 



bo(m) = 1 + — - 5 



IT 



(13) 



Eq. ( [Tl| ) gives subleading effects within an strict expansion in as- Therefore, it can be ap- 
proximated to 



d ^ 



CACf 



Di'^ - {as - isis + l)Dfl - + 4Dg} , (14) 



and the solution reads 



Bsiup) = bs{m) + Ai^^ ln(w;^«) + A2as{m) 



W' 

I3o-13Ca/6 _ I 



+ Asasim 



fSo-2CA 



+ A^asim) ln(2^°) 



(15) 



where Po = yC^ — ^TpUf, z 
Eq. (|15D read 



as{m) 



1 

/3o 



and w 



/3o 



The coefficients Ai in 



A, 

Ay- 
A,-- 
A, 



7TCf[3(3o{26Cl + 19CACf - 32Cj) 



- Ca(208C1 + 651CaC/ + 116C|) 



PoiMs + 1) - 3) + Ca(15 - 14s(s + 1)) 



6(/3o - 2Ca)2 
247rCK3/3o - 11Ca)(5Ca + 8Cf 



13Ca(6/3o-13Ca) 



CI(-9Ca + 100C 



f) 



Pii6Po-l^CAm-2CA) 
+(3oCA{-7ACf + Ca{^2 - 13s{s + 1))) + 6/3^{2Cf + Ca{-3 + s{s + 1))) 



(16) 



Our evaluation can be compared with the result obtained using the vNRQCD formalism 
1^. We agree for the spin- dependent terms but differ for the spin-independent ones. The 
disagreement still holds if we consider QED with light fermions {Cj 1, Ca — ^ 0, Tp — 1). 
Agreement is found if we consider QED without light fermions {Cf 1, Ca 0, Uf ^ 0, 
Tp 1). If we expand our results in as, we can compare with earlier results in the literature. 
By following the discussion in Ref. M, we can relate our results with the correction to the 



wave-function at the origin as defined in Ref. |16|. We obtain 
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BJm] 



2--s(s + l) 



r41 7 



- -sis + 1) CfCA + -C 



12 12 



2 



(17) 



2 



(2--s(s + 1))C/ + Ca 



+ . . . 



where we have expanded up to second order in ln(z/p) = In(mas) with as = asii'p). The first 
term reproduces the leading log term |1T^JT9[| (see also 0]), the /3o-independent 0{a^\n^ a^) 
terms reproduce Kniehl and Penin results |TB| and we agree with the complete 0{alln^ a^) 
term computed by Manohar and Stewart [0] (the sign of difference for the /5o-dependent terms 
displayed in Ref. 0] is due to the fact that in Ref. [0] the expansion was made with as(m) 
whereas here we have chosen as{mas)). Nevertheless, disagreement with this last evaluation 
appears at higher orders in the expansion in as (we have explicitly checked this for the 
0(ag In^tts) terms). As far as we can see, the disagreement seems to be due to the fact that 
they have different expressions for the RG improved potentials PJ^l-^ 



By setting 



Bs{vp) includes all the large logs at NLL order in any (inclusive enough) 



S-wave heavy-quarkonium production observable we can think of. For instance, the decays 
to e^e~ and to two photons at NLL order read 



T{VQ{nS) 



e e 



r(PQ(n5)^77) = 6 
~ 6 





2 


fniQCfas 


.MvQ{nS)_ 




K n 




2 


(mqCfas 


.MvQ{nS)_ 




K n 


C^evaQ 


2 


/niQCfas 


MpQ{nS)_ 




\ n 




2 




MpQ{nS)_ 




\ n 



{Bi{Up){l + 5<Pn)Y 



{l + 2(5i(z/p)-l) + 250„,} , 



{B,{Vp){l + 5<Pn)f 



{l + 2(5o(z/p)-l) + 250„} 



(19) 



where V and P stand for the vector and pseudoscalar heavy quarkonium, we have fixed 
Vp = mgCfas/n, as = as{h'p), and ("ifniz) = ^^^J'^^ and T{z) is the Euler F-function) 



TT 



Ca + ^ f ^i(n + 1) - 2n^2(n) + ^ + 7e + -) 
4 V In) 



(20) 



which has been read from Ref. |]TP[ (see also |2^). Working along similar lines one could 



' The running of the Coulomb potential is not needed for the precision of the above calculation. 



Nevertheless, the running obtained in vNRQCD |21| also disagrees with the one obtained in pN 



RQCD At this respect, we would like to report on a recent computation |22| of the 4-loop 
double log term of the Coulomb potential proportional to C^/3o that agrees with the pNRQCD 
result and disagrees with the vNRQCD one. 
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easily obtain NLL expressions for other heavy quarkonium observables in the study of t-t 
production near threshold or in sum rules of bottomonium. Note that for t-t production near 
threshold there already exists a (partial) NNLL RG improved evaluation within the vNRQCD 
formalism |2^. Since we disagree for the RG improved expression for the electromagnetic 



current matching coefficient, this discrepancy would also propagate to that evaluation. 

In conclusion, by using the method of Ref. p| and incorporating the idea |^ of correlating 
the cut-offs of the effective theory, we have given the first steps towards the creation of 
a comprehensive system of RG equations in pNRQCD once the scale Up enters into the 
game. We have used this formalism to compute the running of the matching coefficients 
of the vector and pseudoscalar currents and disagreement with the results obtained using 
the vNRQCD framework |^ has been found. Our results allow to obtain S-wave heavy 
quarkonium production observables with NLL accuracy. We have explicitly illustrated this 
point for heavy- quarkonium decays to e~^e~ and to two photons. We have also computed the 
0{ma^liii^ a) corrections to the Hydrogen spectrum in the Appendix. 
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1 Appendix: 0{ma^\n^a) contributions to the Hydrogen energy 



With the above discussion, we may also try to see whether we are able to obtain the 
0{ma^\n^ a) contributions to the Hydrogen spectrum. It goes beyond the scope of this 
paper to perform a detailed analysis. Here, we will just see that under some assumptions 
from Ref. [O], we are able to obtain the 0(ma^ln^ a) correction to the Hydrogen spectrum. 



We will use the notation of Ref. ||25|. In that paper. Hydrogen-like atoms with light fermions 



were considered. Here, we will use their results on the strict Hydrogen limit (no light fermions: 
rif = 0). 



According to Ref. [|T^], the 0{ma In a) correction to the Hydrogen energy can be obtained 

in 



from the anomalous dimension due to diagrams of the type of Fig 2. a in Ref. fl^ 
our case, to the diagram in Fig. 1. The argument that led to this conclusion was that the 
0(mQ;® In'^ a) terms had the highest possible log power that could appear from a NNNLL 
evaluation of the energy and that, in order to achieve such power, it was necessary to mix 

(2) 

with NNLL logs. The latter only appear in the LL evaluation of , which, indeed, only 



produces a single log (see also [Q). The other point was that the NLL evaluation of 
the potentials would only produce single logs unless mixed with LL running. Therefore, the 
diagrams with the highest possible power of Df^ will give the highest possible log power of 
the Hydrogen energy at NNNLL. 
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The RG equation for the coefficient of the delta potential due to Fig. 1 reads 



^ n(2)/ \ 



Z^aD 



(2)2, 



(21) 



The loop integral of the diagram of Fig. 1 is just equal to the one that gives the running of 
Cs due to Dd^s- Therefore, we can obtain Eq. (|2lD from Eq. ([TT|) by just introducing a factor 
4 due to the fact that we have to change the reduced mass, m/2, from the equal mass case 
to the reduced mass, m, from the Hydrogen-like case. The left-hand side of Eq. (^) gives 

(2) 

the relevant running of D'^^' with NIL accuracy for our case. Therefore, in the right-hand 



side we only need with LL accuracy, which we read from Ref. 25 1 in the limit nj = 0: 



a 



CD{yl/m) 



(22) 



where 



16 a, z/„ 
m — 

3 TT m 



(23) 



Therefore, Eq. (pT]) approximates to 

The above equation gives the following correction at 0(q;^ In'^): 



(24) 



64 



27 \n) m 



(25) 



This contribution gives the following correction to the Hydrogen spectrum {vp ~ mZa): 



SE 



64 
27 



m{Za) 



a 



TT 



W Za 



(26) 



This result agrees with the analytical result of Karshenboim [EB] , the numerical computation 



of Goidenko et al. and the analytical result of Manohar and Stewart |[T^. It disagrees 

^, which 



with the numerical computations of Malampalli and Sapirstein p8| and Yerokhin 
agree with each other. Nevertheless, it may happen that the latter computations are not 



complete for the desired accuracy [30 . 



References 



[1] W.E. Caswell and G.P. Lepage, Phys. Lett. B167, 437 (1986); G.T. Bodwin, E. Braaten and 
G.P. Lepage, Phys. Rev. D51, 1125 (1995); (E) ibid. D55, 5853 (1997). 

[2] A. Pineda and J. Soto, Nucl. Phys. B (Proc. Suppl.) 64, 428 (1998) 



11 



[3] A. Pineda and J. Soto, Phys. Lett. B420, 391 (1998); Phys. Rev. D59, 016005 (1999). 

[4] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. Rev. D60 (1999) 091502; Nucl. Phys. 
B566, 275 (2000); Phys. Lett. B470 (1999) 215. 

[5] M.E. Luke, A.V. Manohar and I.Z. Rothstein, Phys. Rev. D61, 074025 (2000). 

[6] A.V. Manohar and I.W. Stewart, Phys. Rev. D62, 014033 (2000); Phys. Rev. D62, 074015 
(2000). 

[7] A.V. Manohar and I.W. Stewart, Phys. Rev. D63, 054004 (2001). 
[8] A. Pineda, Phys. Rev. D65, 074007 (2002). 

[9] A.V. Manohar, Phys. Rev. D56, 230 (1997); A. Pineda and J. Soto, Phys. Rev. D58, 114011 
(1998). 

[10] E. Eichten and B. Hill, Phys. Lett. B243, 427 (1990); A.F. Falk, B. Grinstein and M.E. Luke, 
Nucl. Phys. B357, 185 (1991); B. Blok, J.G. Korner, D. Pirjol and J.C. Rojas, Nucl. Phys. 
B496, 358 (1997). C. Bauer and A.V. Manohar, Phys. Rev. D57, 337 (1998). 

[11] A. Pineda and J. Soto, Phys. Lett. B495, 323 (2000). 

[12] A.V. Manohar and I.W. Stewart, Phys. Rev. Lett. 85, 2248 (2000). 

[13] M. Beneke and V.A. Smirnov, Nucl. Phys. B522, 321 (1998). 

[14] M. Beneke, A. Signer and V. Smirnov, Phys. Rev. Lett. 80 2535 (1998); A. Czarnecki and K. 
Melnikov, Phys. Rev. Lett. 80, 2531 (1998). 

[15] R. Barbieri, R. Gatto, R. Kogerler and Z. Kunszt, Phys. Lett. B57, 455 (1975). 

[16] B.A. Kniehl and A. A. Renin, Nucl. Phys. B577, 197 (2000); Nucl. Phys. B563, 200 (1999). 

[17] W. Fischler, Nucl. Phys. B129, 157 (1977); S.N. Gupta, S.F. Radford and W.W. Repko, Phys. 
Rev. D26 3305 (1982); S. Titard and F.J. Yndurain, Phys. Rev. D49, 6007 (1994); Y. Schroder, 
Phys. Lett. B447, 321 (1999); B.A. Kniehl, A.A. Renin, V.A. Smirnov and M. Steinhauser, 
Phys. Rev. D65, 091503 (2002); M. Peter, Phys. Rev. Lett. 78, 602 (1997). 

[18] I. Harris and L.M. Brown, Phys. Rev. 105, 1656 (1957); R. Barbieri, E. d'Emilio, G. Curci and 
E. Remiddi, Nucl. Phys. B154, 535 (1979). 

[19] K. Melnikov and A. Yelkhovsky, Phys. Rev. D59, 114009 (1999); A.A. Renin and A.A. 
Pivovarov, Nucl. Phys. B549, 217 (1999). 

[20] A.H. Hoang and T. Teubner, Phys. Rev. D58, 114023 (1998). 

[21] A.H. Hoang, A.V. Manohar and I.W. Stewart, Phys. Rev. D64, 014033 (2001). 

[22] A.A. Renin and M. Steinhauser, private communication. 

[23] J.H. Kiihn, A.A. Renin and A.A. Pivovarov, Nucl. Phys. B534, 356 (1998); F.J. Yndurain, 
|hep-ph/0007333 . 

[24] A.H. Hoang, A.V. Manohar, I.W. Stewart and T. Teubner, Phys. Rev. Lett. 86, 1951 (2001). 



12 



[25] A. Pineda, |hep-ph/02042T^ , to be published in Phys. Rev. A. 

[26] S.G. Karshenboim, Sov. Phys. JETP 76, 541 (1993). 

[27] I. Goidenko et al, Phys. Rev. Lett. 83, 2312 (1999). 

[28] S. Mallampahi and J. Sapirstein, Phys. Rev. Lett. 80, 5297 (1998). 

[29] V.A. Yerokhin, Phys. Rev. A62, 012508 (2000). 

[30] V.A. Yerokhin, Phys. Rev. Lett. 86, 1990 (2001). 



13 




Fig. 1. One possible contribution to the running of D^ l. at NLL. The first picture represents the 

calculation in terms of the free quark- antiquark propagator g'^^ and the small rectangles the poten- 
tials. The picture below is the representation within a more standard diagrammatic interpretation 
in terms of quarks and antiquarks. The delta potentials are displayed as local interactions and the 
Coulomb potential as an extended in space (but not in time) object. 
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Fig. 2. Diagrams, up to permutations, that contribute to the running of B. 
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